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10 CROSS REFERENCE TO RELATED APPLICATIONS 

This U.S. Patent Application is a continuation of U.S. Patent Application No. 09/790,591 
filed February 23, 2001 entitled "Continuous Measurement- While- Drilling Surveying". 

15 FIELD OF THE INVENTION 

The present invention relates to continuous measurement- while-drilling surveying methods 
and apparatuses. 

20 BACKGROUND OF THE INVENTION 

To obtain hydrocarbons such as oil and gas, boreholes are drilled by rotating a drill bit 
attached at a drill string end. A large proportion of the current drilling activity involves 
directional drilling, i.e., drilling deviated and horizontal boreholes to increase the hydrocarbon 

25 production and/or to withdraw additional hydrocarbons from the earths formations. Modern 
directional drilling systems generally employ a drill pipe having a drill bit at the bottom that is 
rotated by a drill motor or a "mud motor". Pressurized drilling fluid or "mud" or "drilling mud" is 
pumped into the drill pipe to rotate the drill motor and to provide lubrication to various members 
of the drill string including the drill bit. The drill bit and drill motor form part of what is known 

30 as the bottom hole assembly ("BHA"). As required the drill pipe is rotated by a prime mover, 
such as a motor, to facilitate directional drilling and to drill vertical boreholes. 



Measurement- While-Drilling (MWD) surveying for directional and horizontal drilling 
processes is performed to provide the orientation and the position of the BHA [Conti, 1999]. 
Azimuth, the inclination and the tool face angles determine the orientation of the BHA, while 
5 latitude, longitude and altitude determine the position of the BHA. The altitude directly 

determines the true vertical depth of the BHA. State of the art MWD surveying techniques are 
based on magnetic surveying which incorporates three-axis magnetometers and three-axis 
accelerometers arranged in three-mutually orthogonal directions. The three-axis accelerometers 
monitor the Earth gravity field to provide the inclination and the tool face angles. This 
10 information is combined with the magnetometer measurements of the Earth magnetic field to 
provide the azimuth [Figure 1, Russel et al., 1979]. 

The magnetic surveying system determines the BHA orientation at certain survey stations 
with the assumption that the error which modifies the Earth's magnetic field vector at the 
15 surveying instrument is in the direction of the borehole [Russel et al., 1979]. This assumption is 
justified by installing these magnetometers inside a non-magnetic housing. Such housing system 
necessitates the use of non-magnetic drill collars around the surveying equipment at a cost 
approaching $30,000 per single installation [Rehm et al., 1989]. 

20 Although simple, the magnetic surveying system suffers from several inadequacies 

especially within the drilling environment [Thorogood, 1990; Thorogood et al., 1986]. The 
presence of downhole ore deposits deviate the measurements of the Earth magnetic field even 
with the non-magnetic drill collars surrounding the surveying instruments. In addition, magnetic 
surveying tools located in non-magnetic drill collars are subject to the influence of the other steel 

25 components of the drill string. It has been shown that drill string induced surveying error 

increases with inclination and as the borehole direction approaches the east-west azimuth. Drill 
string magnetic interference is particularly noticeable when inclination exceeds 30°. Although it 
has been reported that the effect of the drill string magnetic interference could be reduced 
(mitigated but never entirely eliminated) by running long lengths of non-magnetic materials 



above and below the survey instruments [Grindord et al., 1983], this solution could affect the 
cost benefits of the horizontal drilling technology. The effect of drill string magnetic interference 
as well as the presence of downhole ore deposits can neither be quantified nor compensated. In 
addition to these two effects, geomagnetic influences play an important role in the accuracy of 
5 the magnetic surveying system [Parkinson, 1983]. Geomagnetic influences are defined by the 
variation of the dip angle, the declination and the total magnetic field strength with respect to 
time. The dip angle is the angle between the direction of the Earth's magnetic field and the 
horizontal plane. The declination is the angle between the magnetic north and the true north. It 
was recorded that during any given day at a random location the standard deviations of the dip 
1 0 angle, the declination and the magnetic field strength were 0.3°, 0.9° and 0.3 \x -Tesla 

respectively [Parkinson, 1983]. The variation in the geomagnetic field is quite significant in 
relation to the performance capabilities of the magnetic surveying tools currently used. Therefore, 
geomagnetic effects must be taken into account when considering absolute survey accuracy and 
impose definite limitations on the accuracy levels that can ultimately be achieved. 

15 

Several investigations have been reported to improve the magnetic surveying accuracy. 
Shiells and Kerridge (2000) introduced the interpolated in-field referencing method in which 
absolute local geomagnetic field data is determined from spot measurements of the Earth's 
magnetic field. These measurements are taken at local site which is sufficiently close to the 

20 drilling site so that the measurement data is indicative of the Earth's magnetic field at the drilling 
site but is sufficiently remote from the drilling site so that the measurement data is unaffected by 
magnetic interference from the drilling site. This technique compensates for the geomagnetic 
influences but the measurements taken by the magnetometers still suffer from the drill string 
magnetic interference. In addition, measurements of the geomagnetic field in a local site might 

25 not be applicable in an area that has many drilling rigs. Moreover, such technique adds more 
complexity that may give rise to additional cost. 



Hartmann (1998) introduced a method to improve the accuracy of borehole surveying. 
This method is based on determining the uncertainty in the magnetic field measurements by 
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comparing the measured magnetic field with a corresponding known value at specific location. 
However, this method does not take into account two phenomena. The first phenomenon is the 
increase of the drill string magnetic interference at high inclinations. The second is the presence 
of ore deposits that can be found randomly downhole. These two phenomena add additional 
5 errors that are not taken into consideration when referencing the theoretical value of the magnetic 
field to the measured value at the known location. 

Dipersio (1995) introduced a new method for the compensation of geomagnetic 
influences. This method depends on matching both the calculated magnetic field strength and the 
10 calculated dip angle to their corresponding nominal values at a particular geomagnetic location. 
This results in determining an error-free value of the axial component of the magnetic field 
which is directly used to determine the azimuth. However, at high inclinations, the drill string 
magnetic interference generates additional errors in the cross-axial directions. 

1 5 Several other techniques have been described for magnetic MWD surveying [Nicholson, 

1995; Engebretson, 1992; Helm, 1991; van Dongen et al., 1987; Trowsdale, 1982]. However, it 
is clear that the major weaknesses of the present directional sensing instruments stem from the 
use of magnetometers for monitoring the azimuth and from the hostile environment in which 
these devices work. The problem encountered with the use of magnetometers is the presence of 

20 massive amount of steel around the drilling rig. The abundance of ferromagnetic material 

necessitates the separation of the magnetometers by non-magnetic drill collars. Aside from the 
cost of utilizing non-magnetic drill collars, their use introduces a second problem. Since the non- 
magnetic drill collars impose an additional weight on the drill bit, the surveying tools are 
separated from the bearing assembly and the drill bit by about 30 feet [Conti, 1989; Rehm at al., 

25 1989]. Elimination of the non-magnetic drill collars could reduce the distance between the 

instrument package and the drill bit by approximately 50%. The third problem associated with 
the use of magnetometers is their lack of reliability when used underground due to the deviation 
of the Earth's magnetic field from ore deposits. 
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It is known to use strapdown inertial navigation systems for borehole surveying using 
probes. Borehole surveying is performed after finishing the drilling process, or sections of the 
drilling process, for the purpose of quality assurance. U.S. Patent No. 4,542,647 describes a 
borehole inertial guidance system which utilizes two ring laser gyroscopes having their sensitive 
5 axes normal to the borehole axis. In this case, azimuth is determined by a synthetic rotation 
signal provided by three-axis accelerometers. While this configuration may be suitable for 
borehole mapping, it fails to be useful for MWD surveying. There is considerable rotation about 
the tool-spin axis in a MWD process making the use of a synthetic signal to monitor this rotation 
and therefore azimuth useless. 

10 

In addition, ring laser gyros are unsuitable for MWD processes because they cannot 
withstand the harsh conditions encountered near the BHA during the drilling process. In 
particular, the vibration and shock forces associated with the progress of the drill bit impairs the 
effectiveness of ring laser gyros. 

15 

Therefore, there is a need in the art for a continuous measurement-while-drilling surveying 
method and apparatus which mitigates the difficulties of the prior art. 

SUMMARY OF THE INVENTION 

20 

The present invention relates to the replacement of the MWD magnetic surveying system 
with the technology of inertial navigation. Inertial navigation systems (INS) determine the 
position, velocity and orientation of a moving body in three-dimensional space by integrating the 
measured components of the acceleration and the angular velocity of the body [Titteron et al, 
25 1997]. The measurements of acceleration and angular velocity are provided by three-axis 
accelerometers and three-axis gyroscopes arranged in three mutually orthogonal directions. 

Currently available 3-gyroscopes/3-accelerometers INS technology cannot provide high 
accuracy and small dimensions within the same setup and therefore at present, the direct 
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installation of such system inside the BHA to replace the set of magnetometers is not feasible. 
Thus, the traditional high accuracy INS setup is not applicable downhole due to the limited space 
inside the BHA and the relatively large dimensions of the gyroscopes. 



5 Therefore, herein disclosed is a method and apparatus for downhole strapdown INS 

utilizing only two novel gyroscopes. One gyroscope has its sensitive axis along the tool spin axis 
while the other has its sensitive axis normal to the tool spin axis and along the forward direction. 
This results in the removal of the costly non-magnetic drill collars because the gyroscopes are 
insensitive to magnetic fields. In one embodiment, the two gyroscopes are mounted inside the 

10 bearing assembly and are fiber optic gyroscopes (FOG). FOGs are suitable for the downhole 

drilling application because they do not contain moving parts, thus providing high reliability with 
no need for frequent calibration or maintenance [Kim, 1994; Merhav, 1993; Lefevre, 1993]. In 
addition, the FOGs exhibit low environmental sensitivity since they can withstand relatively high 
temperatures, shocks and vibrations [Noureldin et al., 1999; Noureldin et al., 2000]. Moreover, 

1 5 currently available FOGs are of small size, with a drift rate less than l°/hr, long mean time 
between failure (60,000 hrs), no gravitational effects and excellent immunity to vibration and 
shock forces. 

In one embodiment, the FOG that has its sensitive axis along the tool spin axis is 
20 prepared as torus to accommodate the flow of drilling fluid through the bearing assembly. The 3 
accelerometers - 2 FOG system is capable of providing continuous MWD surveying during the 
whole drilling process, up until the horizontal section of the well, instead of MWD station-based 
surveying. In addition, continuous MWD surveying is possible through an initial radical section 
of the well with 3 accelerometers and a single FOG which has its sensitive axis aligned with the 
25 tool spin axis. The initial radical section of the well may be that portion with an inclination less 
than about 45°, preferably less than about 30° and more preferably less than about 20°. 
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Therefore, in one aspect, the invention comprises a continuous measurement-while-drilling 
surveying apparatus for surveying the drilling progress of a bottom hole assembly ("BHA") 
comprising: 

5 (a) a housing adapted to be part of the BHA, said housing having a tool-spin axis; 

(b) a first fiber-optic gyroscope mounted within the housing for generating a first angular 
rotation signal representing angular rotation of the BHA about the tool-spin axis; 

(c) a second fiber-optic gyroscope mounted within the housing for generating a second 
angular rotation signal representing angular rotation of the BHA about an axis normal 

10 to the tool-spin axis; 

(d) accelerometer means for generating three acceleration signals representing the 
components of acceleration of the BHA along three mutually orthogonal axes; 

(e) first processing means responsive to the acceleration signals for determining the 
angle of the BHA away from the vertical and for generating a third angular rotation 

1 5 signal representing rotation of the BHA about an axis normal to the sensitive axes of 

the first and second gyroscopes; 

(f) second processing means responsive to the first, second and third angular rotation 
signals and the acceleration signals for transforming signals representing movement 
of the BHA in a BHA coordinate system to a earth local-level coordinate system; 

20 (g) third processing means operatively connected to the second processing means for 

determining the orientation of the BHA, determining the velocity changes of the 
BHA, updating the velocity components of the BHA and updating the position 
components of the BHA. 

25 The first gyroscope may be toroidal. The first, second and third processing means may be a general 
purpose computer programmed with appropriate software, firmware, a microcontroller, a 
microprocessor or a plurality of microprocessors, a digital signal processor or other hardware or 
combination of hardware and software known to those skilled in the art. 
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In another aspect of the invention, the invention is a method of continuous MWD surveying 
of a wellbore which includes, on completion, a vertical section, a radical section and a horizontal 
section, using a BHA comprising a first gyroscope having its sensitive axis aligned with the tool spin 
axis, a second gyroscope having its sensitive axis on an axis normal to the tool spin axis and three 
5 accelerometers arranged in three mutually orthogonal directions, said method comprising the steps 
of: 



(a) receiving the rotation signals from the first gyroscope and the second gyroscope and 
the acceleration signals from the three accelerometers; 
10 (b) establishing the desired azimuth direction; 

(c) determining the pitch angle from the accelerometer signals; 

(d) determining the rotation rate about an axis normal to the sensitive axes of the two 
gyroscopes by determining the time rate of change of the pitch angle; 

(e) determining the effect of the earth rotation and local-level (/-frame) change of 
1 5 orientation to the angular changes monitored along the three axes and removing these 

two effects from the rotation rate signals; 

(f) transforming the rotation signals from a BHA coordinate frame to the earth local- 
level coordinate frame; 

(g) determining the pitch, roll and azimuth of the BHA; 

20 (h) determining the time rate of change of the velocity of the BHA expressed in the £ - 

frame and update the velocity components of the BHA; and 
(i) updating the position components of the BHA. 



In one embodiment, this method is used for an initial portion of the radical section of the 
25 wellbore. The initial portion of the radical section of the well may be that portion with an 
inclination less than about 45°, preferably less than about 30° and more preferably less than 
about 20°. After the initial portion, the method is altered by switching the axes orientation of the 
two gyroscopes. Once the wellbore reaches the horizontal section and rotary mode drilling is 
required, the surveying method may then switch to station-based surveying. 
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BRIEF DESCRIPTION OF THE DRAWINGS 



Figure 1 is a block diagram illustrating the functionality of a prior art magnetic surveying 
5 system. 

Figure 2 is a sketch representing the inertial frame. 

Figure 3 is a sketch representing the Earth-fixed frame. 

10 

Figure 4 is a sketch representing the local level frame. 

Figure 5 is a sketch showing the body frame axes arranged inside the bearing assembly of a 
BHA. 

15 

Figure 6 shows the change of t -frame orientation along the Earth's surface. 

Figure 7 shows the arrangement of the fiber optic gyroscopes inside a bearing assembly. 

20 Figure 8 is a bottom view for the bearing assembly showing the flow of mud through the 

toroidal FOG. 

Figure 9 shows the initial arrangement of the FOGs and the accelerometers with the body 
frame axes inside the bearing assembly. 

25 

Figure 10 shows the setup while surveying the radical section until 20° inclination. 

Figure 1 1 shows the setup after achieving 20° inclination. 

30 Figure 12 shows the new axis orientation after 20° inclination. 

Figure 13 shows the relationship between the pitch and the inclination angles before and 
after changing the axes orientation. 

3 5 Figure 1 4 is a block diagram describing the functionality of the surveying procedure until 20° 

inclination. 

Figure 14A is a block diagram describing the functionality of the surveying procedure 
until 20° inclination using only the toroidal FOG. 

40 

Figure 1 5 is a block diagram describing the functionality of the surveying procedure after 20° 
inclination. 
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Figure 16 shows the arrangement of the surveying instruments inside the bottom hole 
assembly while drilling the horizontal section of the well. 

Figure 17 is a block diagram showing the station-based surveying procedure. 

5 

Figure 18 is a block diagram showing the technique of aided inertial navigation. 
DESCRIPTION OF THE INVENTION 

10 The present invention provides for a continuous measurement-while-drilling surveying 

method and apparatus. All terms not defined herein have their common art-recognized meanings. 

Strapdown Intertial Navigation 

15 There are two common INS categories [Titteron et al., 1997; Britting, 1971]. The first 

category is the gimbaled INS that maintains a stationary platform for the accelerometers by 
torque motors. These systems use the angular rate measurements of the gyroscopes attached to 
the platform to control the motors. The second category is the strapdown INS (SENS) that 
mathematically transforms the output of the accelerometers attached to the body into a locally- 

20 leveled coordinate system before performing integration. These systems use the output of the 
gyroscopes attached to the body to continuously convert from body coordinate to locally leveled 
coordinate. 

Three reference frames are commonly used in the technology of inertial navigation 
25 [Schwarz, 1999; Salychev, 1998]. These are the inertial frame, the earth- fixed frame and the 

local level frame. The inertial frame (i-frame) is a non-rotating and non-accelerating frame. The 
origin of the inertial frame is at the center of mass of the Earth, its Z-axis is parallel to the spin 
axis of the Earth, its X-axis points toward the mean equinoctial colure and its Y-axis completes a 
right handed orthogonal frame (see Figure 2). The Earth-fixed frame (e-frame) has its origin at 
30 the center of mass of the Earth with its Z-axis parallel to the spin axis of the Earth, X-axis 
pointing towards the mean meridian of Greenwich and Y-axis completing a right handed 
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orthogonal frame (see Figure 3). The local level frame (£ -frame) has its origin at the origin of the 
sensor frame with its Z-axis orthogonal to the reference ellipsoid of the Earth and pointing 
upward, Y-axis pointing toward the geodetic north and X-axis completing a right handed 
orthogonal frame (see Figure 4). When the inertial sensors (gyroscopes and accelerometers) are 
5 installed inside a moving object, they are oriented in three mutually orthogonal directions known 
as the body frame (b-frame). The axes of the b-frame point toward the forward (Y), transverse 
(X) and vertical (Z) directions (see Figure 5). 

The measurements of linear acceleration and angular velocities are taken in the b-frame. 
10 These measurements should be transformed into the corresponding reference frame within which 
they are processed to provide the position vector, the velocity vector and the attitude components. 

The transformation from the b-frame to the £ -frame is performed with the rotation matrix r[ , 
which is expressed as a function of the attitude components (azimuth v|/ , roll (tool face) <|> and 
pitch (inclination) 0 ) as follows: 



15 



f cos y/ cos <p - sin y/ sin G sin (j> - sin y/ cos 6 cos y/ sin <f> + sin y/ sin 0 cos ^ 
sin y/ cos (/> + cos y/ sin 6 sin (f> cos y/ cos 0 sin y/ sin <f> - cos y/ sin 6 cos ^ 



-cosflsin^ sin# cos#cos^ 



(1) 



J 



The position of the BHA as described in the I -frame is expressed in terms of the 
curvilinear coordinates (latitude cp , longitude X and altitude h ). The velocity is alternatively 

20 expressed by three components along the East direction ( V e ), north direction ( V n ) and vertical 

direction (V u ). Therefore, the time rate of change of the position components is related to the 
velocity components as follows [Schwarz et al., 1990]: 



-11- 









0 

1 


1 


0 






(<p) 




M + h 




X 

w 




0 
0 


0 

1 


F" 

v / 






{N + h)cos<p 
0 



(2) 



where 

N : is the prime vertical radius of curvature of the Earth's surface (East- West); 
M : is the corresponding meridian radius of curvature (North-South); 



D 



_! : is a 3-by-3 matrix defined as 
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The acceleration of the bearing assembly is measured in three mutually orthogonal 
directions in the b-frame f b =(f x f y f z ) T . These measurements are transformed into the £ - 

5 frame using the rotation matrix (f e = Rlf b ). However, the acceleration components 

expressed in the £ -frame, /* , can not directly give the velocity components of the bearing 
assembly due to three reasons. The first reason is due to the rotation rate of the Earth about its 
spin axis ( co e =1 57hr) which is interpreted in the £ -frame as co\ e = (o co e cos cp a> e sin cp) T . 

The second reason is the change of orientation of the £ -frame with respect to the Earth as the 
10 bearing assembly penetrates the downhole formation. The £ -frame change of orientation along 
one of the meridians can be shown on Figure 6. This change of orientation is due to the definition 
of the local North and vertical directions. The North direction (N) is tangent to the meridian all 
the time while the vertical direction (V le ) is normal to the Earth's surface. This effect is 

interpreted as rotation rate vector expressed as a>^ = (- <p X cos q> X sin <p) T which can be 



15 rewritten as a> et 



-V V e ^ c tan^ 
M + h N + h N + h 



T (see equation 2). The third reason is the 



12 



gravity field of the Earth which is given as g l = (0 0 - y) T , where y is obtained from the 
well-known normal gravity model given as follows [Schwarz et al., 1999]: 

y = a x (1 + a 2 sin 2 <p + a 3 sin 4 + 0 4 + a 5 sin 2 cp)h + <? 6 /z 2 (3) 

where 

ai = 9.7803267715 m/sec 2 34= -0.000003087691089 1/sec 2 

a 2 = 0.0052790414 a 5 = 0.000000004397731 1/sec 2 

a 3 = 0.0000232718 ao = 0.000000000000721 l/(msec 2 ) 

Then, the time rate of change of the velocity components V of the bearing assembly 
while penetrating downhole can be expressed as follows [Schwarz et al., 1999]: 



K £ =^/*-(2n£ + n^+^ (4) 

where Qf e and Q £ ei are the skew-symmetric matrices corresponding to co\ e and co E ei respectively 
and they are expressed as follows: 
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(6) 



The rotation matrix R b can be obtained by solving the following differential equations 
[Schwarz et al., 1999]: 
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(7) 



where Q] b is the skew-symmetric matrix of the measurements of angular velocities provided by 

the gyroscopes. Q b ib corresponds to the angular velocity vector co\ h = (co x co y a> 2 Y and is 
given as: 



4 = 



0 

a, 



2 OJ y 

0 -0)y 



(8) 



Since the gyroscopes measure both the Earth rotation and the change in orientation of 
the t -frame in addition to the angular velocities of the bearing assembly, the term Q b ( in 

h h 

10 equation 7 is subtracted from Q° b to remove these two effects. The term Q i( consists of two 

parts. The first part is which accounts for the Earth rotation rate and the second part is Q. b e 

which accounts for the orientation change of the I -frame. Therefore, Cl b it can be written as 
follows: 



15 



(9) 



Since a h ie = R b co\ e and co b ( = R b (o l eSL , then co b e = R b {co l ie + co l el ) and can be given as 



follows: 
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Consequently, the corresponding skew-symmetric matrix Q h it can be determined. 



When equations 2, 4 and 9 are combined together, they form what is known as the 
5 mechanization equations in the I -frame and they are usually given together as follows [Schwarz 
et al., 1990; Schwarz et aL, 1999]: 
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10 The inputs to these mechanization equations are the gyroscope and accelerometer 

measurements while the outputs are the curvilinear coordinates, three velocity components and 
three attitude components. The £ -frame is chosen to describe the drill bit motion for three 
reasons. The first is that the axes of the £ -frame are aligned towards the local East, North and 
vertical (normal to the Earth's surface and pointing upward). Using this orientation, the pitch, the 

15 roll and the azimuth angles describing the drill bit attitude (inclination, tool face and azimuth) 
can be obtained directly as outputs of the mechanization equations solved in the I -frame. The 
second reason for choosing the t -frame is that the geodetic coordinates expressing the drill bit 
latitude, longitude and depth can be obtained directly by solving the mechanization equations 
described in this frame. Finally, the computational errors in the navigation parameters on the 

20 horizontal plane (North-East plane) are bounded due to Schuller' s effect. This effect stipulates 
that the inertial system errors of the horizontal plane components are coupled together producing 
what is called a Schuller loop [Mohammed, 1999; Titterton et al., 1997]. Therefore these errors 
oscillate with frequency called the Schuller frequency (1/5000 Hz). 



25 Apparatus 
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At present, the FOG size determines the overall size of the inertial measurement unit. As 
mentioned earlier, present technological level of FOG-based inertial measurement units do not 
permit direct installation of 3-FOGs inside the BHA. Therefore, the present invention comprises 
a novel structure of inertial measurement unit with two specially designed FOGs in order to 
5 accommodate space and precision requirements. These two FOGs monitor two rotations rate 
signals along the Y and Z directions of the b-frame axes. The third rotation rate signal may be 
generated using the time rate of change of the inclination angle calculated with the three-axis 
accelerometer measurements. The two FOGs and the three-axis accelerometers may be mounted 
anywhere inside the drill pipe and the technique provided in the following sections can process 
10 their measurements to determine the surveying data. 

In one embodiment, the MWD inertial surveying sensors are fitted inside the bearing 
assembly close behind the drill bit as is shown in Figure 7. The system comprises a first toroidal 
FOG (10), a second FOG (12), three accelerometers (14, 16, 18), a microprocessor (20), wireless 
communication transmitter (22) and a power source such as a battery. The microprocessor (20) 
receives the output signals from the two FOGs (10, 12) and the three accelerometers (14, 16, 18) 
and determines the inertial navigation components as described herein. These components are 
then transmitted to the surface for review by a drilling engineer. Alternatively, some or all of the 
processing may take place at the surface if the output signals of the inertial sensors are 
transmitted to the surface. Standard transmission or telemetry apparatus and methods may be 
used, as is well known in the art of MWD. 

The toroidal FOG (10) is installed with its sensitive axis along the tool spin axis. It's 
shape allows the flow of mud through the drill pipe as is clearly shown in Figure 8. The 
25 dimensions of the toroidal FOG (10) are dictated by the dimensions of the bearing assembly. 
The second FOG (12) is installed so that its sensitive axis is normal to the tool spin axis and 
towards the forward direction (Y-axis) Since its sensitive axis is normal to the tool spin axis, this 
FOG is denoted as the normal FOG (12). Both FOGs should preferably have drift rate less than 
1 deg . / hr and angle random walk of less than 0. 1 deg . / hr 1 4Hz in 1 00 Hz bandwidth. The 
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15 



20 



bias drift is defined as the deviation in the measured rotation rate at constant temperature. The 
angle random walk is the broadband random noise component at the FOG output that results 
either from the shot noise or the thermal noise in the photodetector. 

5 The three-axis accelerometers (14, 16, 18) are arranged in three mutually orthogonal 

directions with one of them having its sensitive axis parallel to the tool spin axis. 

Standard wireless communication between the system inside the bearing assembly and 
the surface MWD processor is provided [Skillingstad, 2000]. A power source such as a battery 
10 package provides the necessary power supply to the two FOGs, the accelerometers, the 
microprocessor and the wireless communication system. 

Once the vertical section of the well is established, the horizontal drilling process 
involves three main tasks: 

15 

1 . Establishing the desired azimuth direction while the drill pipe is still in the vertical 
direction. 

2. Drilling the radical section of the well using steering mode of operation. 

3. Drilling the horizontal section of the well with using rotary mode of operation. 

20 

The establishment of the desired azimuth direction 

In order to establish accurately the azimuth direction, precise value of the initial azimuth 
should be known. This initial azimuth value is determined when the whole system is completely 
25 stationary. Therefore, the two FOGs monitor the Earth rotation, while the accelerometers monitor 

the Earth gravity. The Earth rotation vector monitored by the gyroscopes, co b ie , and the Earth 
rotation vector expressed in the t -frame, co\ e , are related through the rotation matrix 
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Rt = {rIJ which transforms the ^ -frame to the b-frame. This relationship can be written as 
follows: 



(o£=(a) x co y co z ) T ^R b t (ol 

'cos^cos^-sin^sin^sin^ sin^cos^ + cos^sintfsin^ ~cos0sin^Y 0 
= -siny/cos# cos^cos# sin# 

^cos^sin^ + sin^sin#cos^ sin^sin^-cos^sintfcos^ cos#cos^ } 



o) e cos<p 
{a) 6 sing) J 



(12) 



The body frame axes are aligned as shown on Figure 9 with the normal FOG having its 
sensitive axis along the Y-axis and the toroidal FOG having its sensitive axis along the Z-axis 
(the tool spin axis). Using equation 12, the rotation rate measurement provided by the normal 

FOG, 0) y , can be given as function of the Earth rotation rate co € as follows: 

a>y = (cos y/ cos 8) of cos#> + (sin#) a> e sirup (13) 
Consequently, the azimuth can be determined as follows: 



sm 9 tan <p 

co e cos q> /t .. 
cost? 



Since the bottom hole assembly is almost in the vertical direction, the inclination angle 
(i.e. the pitch angle, 8 ) is very close to zero. Therefore, we can simplify equation 14 by 
considering sinG « 0 and cos8 » 1 . The initial azimuth can then be given as: 



(O COS (p 
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It should be noted that equation 15 has a singularity when the latitude angle cp equals to ±7t/2. 
This situation corresponds to drilling at the Earth poles that can not exist. 



5 Therefore the initial azimuth can be determined using equation 15 by substituting the 

angular velocity output of the normal FOG, co y , after removing the offset and the constant drift 

values. The residual errors in the measurements provided by this FOG, 8co y , affect the 
calculations of cosy and consequently the accuracy of the initial azimuth. The error in 
cosy due to these errors 8co y can be obtained by differentiating equation 15: 

10 

5co 

£cos^ = — (16) 

co e cos <p 



It is crucial to remove the error 5 cosy from the cosy term calculated from equation 15 before 
proceeding to calculate the corresponding value of the azimuth y . This is necessary because of 
1 5 the inversely proportional relationship between the error in the azimuth 5y and the azimuth 
angle itself, y , which can be obtained from equation 16 as follows: 



( 8(0, 



\* cosy) 
sin^ 



(17) 



20 Since siny < 1 , the term (1/sin^) in equation 17 tends to magnify the error 8co y . When 

y is close to ;r/2 , i.e. approaching the East direction, the term l/siny becomes closer to unity 
and has minor effect on 5y . On the other hand, when y is close to 0, i.e. approaching the North 
direction, the term 1/sin^ magnifies the random error 5y . In order to avoid this problem, the 
error term 5 cosy is determined first using equation 16 and removed from the values obtained 
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for cosy using equation 15. Once an accurate estimation of cosy is determined, the azimuth 
can be computed independently from the drilling direction. 



Once the initial azimuth of the bottom hole assembly is determined, the drill pipe can be 
5 rotated to achieve the desired azimuth direction. Although the normal FOG is responsible for the 
determination of the initial value of the azimuth, the toroidal FOG is responsible for the 
establishment of the desired azimuth direction. 



With the BHA oriented along the vertical direction, the drill pipe is rotated along the tool 
10 spin axis until the desired azimuth direction is established. Rotation along the tool spin axis is 
preferably the only motion, either angular or translational, performed during this operation. 
Under these circumstances, the toroidal FOG monitors the rotation rate along the tool spin axis in 

addition to the component of Earth rotation along the vertical direction (co e sirup). Therefore, the 
time rate of change of the azimuth, i|/ , can be written as follows: 

15 

\\i = co z -co e %\xvq> (18) 



where co z is the measurement delivered by the toroidal FOG. 



20 Equation 18 is solved numerically in real-time to provide continuous monitoring of the azimuth 
angle while rotating the drill pipe. This first-order differential equation can be solved numerically 
using the Euler method [Yakowitz et al, 1989] after substituting \ji by - y/ k )/(A/) , where 

\\i £ + i is the azimuth value at time t k+l , v|/ k is the azimuth value at time t k and A t = t k+l -t k . 

Then, the azimuth can be continuously monitored using the following equation: 

25 

M>*+i =Vk + (co z -a) e sirup) At (19) 



-20- 



As the drill pipe continuously rotates about its tool spin axis, the azimuth angle is 
continuously determined from equation 19. Once the desired azimuth direction is established, the 
rotation of the drill pipe is stopped and building the radical section of the well is initiated. 

5 Surveying the radical section of the well 

During the building of the radical section of the well, the bottom hole assembly starts the 
drilling process from complete vertical direction and continues toward a complete horizontal 
direction. At the beginning, while the BHA is still in the vertical direction, the body frame axes 

10 are chosen as shown on Figure 9. The X, Y, and Z axes point toward transverse, forward and 
vertical directions respectively. Therefore, the Z-axis points toward the tool spin axis. With this 
choice of axes orientation, we assume that the motion of the BHA is along the Y-axis, and 
sliding along the Z-axis (Figure 10). Hence, the toroidal FOG monitors the changes in the 
azimuth angle \y while the normal FOG monitors the changes in the roll angle <|> . However, as 

15 the wellbore inclination increases, the Z-axis approaches the horizontal and the Y-axis 

approaches the vertical. The Y-axis is no longer in the forward direction and the Z-axis is no 
longer in the vertical direction. In addition, the toroidal FOG is no longer in the horizontal plane 
to successfully monitor the azimuth. Therefore, at a predetermined inclination such as 20°, a new 
orientation of the axes may be set to keep Y along the forward direction and Z along the vertical 

20 direction. This can be achieved by virtually choosing new axes system X 1 , Y 1 , and Z\ as is 

illustrated in Figure 12. X 1 points toward the transverse direction exactly like X while Y l points 
toward the forward direction (i.e. along the tool spin axis) and Z 1 points toward the vertical 
direction. This axes orientation is maintained until the end of the drilling process without any 
physical change in the mounting of the two FOGs. 

25 

In a preferred embodiment, there are two surveying methodologies for the radical section of 
the well. The first is applied during an initial portion of the radical section. As used herein, an 
"initial portion" is that portion of the well where the inclination is less than about 45°, preferably less 
than about 30° and more preferably less than about 20° inclination. The body frame is defined by the 
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X, Y and Z axes as is shown in Figure 10. During this initial portion of the well, the toroidal FOG 
(with its sensitive axis along the Z-axis) is responsible for monitoring the changes in the azimuth 
angle while the normal FOG (with its sensitive axis along the Y-axis) is responsible for monitoring 
the changes in the roll angle. Since the toroidal FOG monitors the rotation along the tool spin axis, 
5 the tool face angle is directly obtained from the toroidal FOG output. 

After the initial portion, the tasks of the two FOGs are switched. In addition, the roll and the 
tool face angles become the same and are monitored by the toroidal FOG while the normal FOG 
monitors the changes in the azimuth angle. 

10 

It should also be noted that the relationship between the pitch angle and the inclination angle 
is not the same for the two cases during and after the initial portion of the well. Before 20° 
inclination the pitch angle is equal to the inclination angle while after 20° inclination there is njl 
radians phase shift between them. This is because the pitch angle is defined as the angle between the 
15 forward direction (Y-axis) of the moving body and the horizontal plane (Figure 13), Since the 
definition of the Y-axis is changed at 20° inclination, the definition of the pitch angle is changed 
accordingly and the njl radians phase shift is created between the inclination and the pitch angle. 

The pitch angle is determined from the three-axis accelerometers. The determination of the 
20 inclination angle is based on the assumption that the three-axis accelerometers monitor the gravity 
components in three mutually orthogonal directions and their outputs can be integrated together to 
determine the pitch. This assumption is valid as long as the penetration rate through the downhole 
formation is relatively small. The average angle-building rate is about 4°/100 ft. for long-radius 
application, 12°/100 ft. for medium-radius applications and 4°/ft. for short-radius applications. In 
25 each of these three categories of horizontal drilling, the above assumption can be applied and the 
pitch can be determined from the measurements of the accelerometers. 

Surveying the radical section before 20° inclination 
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The surveying technique explained in this sub section provides continuous surveying 
based on the measurements delivered by the three-axis accelerometers, the toroidal FOG and the 
normal FOG. The toroidal FOG determines the rotation rate along the Z-axis, co z , while the 
5 normal FOG determines the rotation rate along the Y-axis, a> y . The rotation rate along the X- 
axis, co x , can be determined from the time rate of change of the pitch angle. Therefore, the 

process starts with receiving the accelerometer measurements that can be combined to determine 
the pitch, which, until 20° inclination, is the same as the inclination angle. Based on the 
assumption of very small penetration rate through the downhole formation, we assumed that the 
10 three accelerometers are basically affected by the gravity field of the Earth only. Therefore, the 

accelerometer measurements taken at the body frame, g b , are related to the gravity field of the 
Earth expressed in the I -frame, g , through the rotation matrix R° as follows: 



15 





(f x ) 


f 




fy 


-K e g - 




UJ 


V 



cos^cos<z) - sin y sin # sin 0 sin^cos^ -f-cos^sin^sin^ - cos#sin^ 

-siny/xostf cosy/xos# sin# 

^cos y/ sin (j> + sin y/ sin 0 cos (f> sin y/ sin (j> - cos y/ sin 9 cos <j> cos 0 cos ^ 



'0^ 
0 



(20) 

where f x , f y and f z are the measurements provided by the three mutually orthogonal 

20 accelerometers. 

y is the value of gravity obtained from the normal gravity model given in equation 3. 



25 



From equation 20, the accelerometer measurements f x , f y and f 2 can be given in terms of y 
as follows: 
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f x = (cos 9 sin <f>) y 
f y = -(sin0) 7 
f z = - (cos# cos^) / 



(21) 
(22) 
(23) 



where 9 is the pitch (inclination) angle and <j> is the roll angle. 

It is clear that the pitch angle can be obtained directly from equation 22 as long as y is 
determined from the normal gravity model (equation 3). The determination of y depends on the 
available information for the latitude cp and the altitude h . However, since the normal gravity 
model is an approximate model, and in order to prevent the propagation of errors, equations 21 
through 23 are manipulated to remove the gravity term and to directly give the pitch angle. From 
equations 21 and 23, we can write 



Equation 25 gives the pitch angle directly in a way that is completely independent of the 
gravity field of the Earth. Then the rotation rate along the X-axis, co x , is determined using the 

time rate of change of the pitch angle as follows: 




(24) 



Then the pitch angle can be obtained by dividing equation 22 by equation 24. 




(25) 



co x = 




(26) 
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where 6 k and 8 k+] are the values of the pitch angle at time t k and t k+] respectively. 



10 



The surveying procedure delivers the navigation data (the position and the orientation) for 
the bottom hole assembly by solving the set of first order differential equations given in equation 

1 1 . The first step in solving these equations is the parameterization of the rotation matrix 

[Schwarz et al., 1999; Salychev, 1998; Titterton et al., 1997]. One method for this purpose is the 
quaternion approach. According to Euler's Theorem [Salychev, 1998], the rotation of a rigid 
body (represented by the b-frame) with respect to a reference frame (represented by the I -frame) 
can be expressed in terms of the rotation angle 0 about a fixed axis and the direction cosine of 
the rotation axis that defines the rotation direction. Thus, quaternion parameters 

Q = (#i 02 03 04 ) T are introduced to describe the rotation of the b-frame with respect to 
the I -frame and they are expressed as follow: 



Q- 





f 










J4J 


V 



(0j0)sin(©/2)^ 

(e y /e)sin(e/2) 

(0 z /0)sin(0/2) 
cos(0/2) 



(27) 



r~2 — 5 — 2 0 0y 0 

1 5 where 0 = J@ x + @ + ® z is the rotation angle and — , — , — are the three direction 

v y - 0 0 0 

cosines of the rotation axis with respect to the I -frame. 



The definition of the quaternion parameters described in equation (27) implies that the 

four quaternion components are not independent because qf + q\ + q\ + q\ = 1 . This means that 

20 only three independent quaternion components are sufficient to describe the rigid body rotation. 
The time rate of change of the quaternion is described by the following first-order differential 
equation: 
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where Q(co) is a skew-symmetric matrix given as follows: 



(28) 



15 



Q(<y) = 



( 0 
-co 2 
co y 

-CD r 



(O z 

0 

-<W, 



6>r 



CO, 



0 



(29) 



5 where o) x ,a> yy co z are the angular velocities of body rotation which are determined by equation 
26 for co x and are monitored by the normal FOG and the toroidal FOG for co y and co z 
respectively, after compensating for the Earth rotation and the I -frame change of orientation. 

To solve the set of first order differential equations given in equation 28, Euler method 
10 can be used to determine the quaternion parameters Q k+l at time tk+i based on the values of the 

quaternion parameters Q k at time t k as follows: 



V2 



A* 



(30) 



where At -t M -t k . 



Once the quaternion parameters are determined at certain time, the rotation matrix can be 
obtained using the following direct relationship: 



20 
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Ri i R 



l 2l 
(31) 
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On the other hand, the quaternion parameters can be obtained from the rotation matrix R b as 
follows: 
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( 0.25(i? 32 








0.25(^,3 




<?3 




0.25(^2i 


-R n )lq A 






{0.5j\ + R u 


+ R 22 + R 3i j 



(32) 



The quaternion parameters are introduced for the parameterization of the rotation matrix 
Rl for three reasons. The first reason is that only four differential equations are solved 

numerically instead of six differential equations if we manipulate the rotation matrix R b directly. 

The second reason is that the quaternion solution avoids the singularity problem that might exist 
10 with some other solution methods. The third reason is the computational simplicity introduced by 
the quaternion. 



One embodiment of the method may be described as follows: 



15 1 . Receive the measurements from the two FOGs and the three accelerometers. Some modern 
FOGs deliver angular increments at their output instead of rotation rate. Similarly, the 
accelerometers deliver velocity increments instead of specific force (acceleration) 
measurements. The angular changes provided by the two FOGs should be compensated for 
the drift rate d as follows: 



20 



A0 V =A0 v -d v At 

y / " (33) 
A0 Z =A0 z -d z At 

where 

A# is the raw measurement of the angular increments provided by the normal FOG; 
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A 9 y is the measurement of the angular increments provided by the normal FOG after compensating 
for the drift rate d y ; 

is the raw measurement of the angular increments provided by the toroidal FOG; 

A <9 Z is the measurement of the angular increments provided by the toroidal FOG after compensating 
for the drift rate d z ; 

At is the difference between the current time sample t k+] and the previous time sample t k . 

Similarly, the velocity increments provided by the accelerometers should be 
compensated for the bias terms b as follows: 

AVj = APj - b x At 

5 Av y =Av y -b y At (34) 

Av z = Av z - b z At 

where 

Av x is the raw measurement of the velocity increments provided by the accelerometer installed 
along the X-axis; 

AVj is the measurement of the velocity increments provided by the accelerometer installed along the 
X-axis after compensating for the bias b x ; 

Av^ is the raw measurement of the velocity increments provided by the accelerometer installed 
along the Y-axis; 

Avy is the measurement of the velocity increments provided by the accelerometer installed along the 
Y-axis after compensating for the bias b y ; 

Av 7 is the raw measurement of the velocity increments provided by the accelerometer installed 
along the Z-axis; 

Av 2 is the measurement of the velocity increments provided by the accelerometer installed along the 

Z-axis after compensating for the bias b z ; 
At is the difference between the current time sample t k+l and the previous time sample t k . 
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2. Determine the pitch angle from the accelerometer measurements as given in equation 25. 
This equation is slightly modified to suit the velocity increment output Av from the 
accelerometers. This modification is based on the direct relationship between the specific 
5 force measurement / , representing the acceleration, and the velocity increments Av given as 
/ = Av/A/ . Thus, the pitch angle, 0 , is determined as follows: 



6 - arctan 



Av 



Av^+Av* 



(35) 



10 Then, the angular changes along X-direction, A 6 X , can be obtained from equation 26 as 

follows 



A0 x =<o x At = [e(t M )-0(t k )] 



(36) 



15 3. Determine the effect of the Earth rotation and the t -frame change of orientation in the 

angular changes monitored along the three axes. The term responsible for the Earth rotation 

and the I -frame change of orientation is a>f e and is given as follows (see equation 10): 



20 



4 =*<('*) 



M + h 

N + h 

V e {t k )\axi<p e . 

— — + (o sine? 

N + h 



(37) 



Then the angular changes corresponding to co^ can be determined as: 
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(38) 



Consequently, the vector of measured angular increments 9% = \A6 X A6 y A6 Z ) T is 

compensated for 0$ to determine the actual bottom hole assembly angular changes d\ h as 
follows: 

e b tb =9t b -et ( ={se x se y 59 2 ) T (39) 

4. Update the quaternion. The time rate of change of the quaternion as given in equation 28 
corresponds to the time rate of change of the rotation matrix r[ as given in equation 7. The 
quaternion update follows equation 30. The time interval A / is multiplied by the elements of 
the skew-symmetric matrix of angular velocities Q , so that co x At = 86 x , co y At- 86 y and 
co 2 At = 86 z . The quaternion update equation is then given as follows: 
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(40) 



where the initial value of the quaternion Q(t 0 ) = (q { (t Q ) # 2 (*g ) 03 ('o ) #4 ('o )) f 1S chosen to 
be equal to (0 0 0 

Consequently, the rotation matrix R b can be determined from its direct relationship to the 
quaternion Q = (q x q 2 q^ q^) T as given by equation 3 1 . From the definition given in 
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equation 1, the pitch (inclination) is determined as 6 = arcsin(# 32 ) , the roll is determined as 
<j) = arctan(- #31/^33) and the azimuth (tool face) is determined as y = arctan(- #12/^22) 

5. Determine the changes in velocity of the bottom hole assembly using equation 4 that 
5 describes the time rate of change of the velocity of the bottom hole assembly expressed in the 
£ -frame. This equation can be rewritten as follows: 

At 

or (41) 
AV i ^Rlf b M-(2Qi +n e ei )V e At + g £ At 

10 Since the specific force measurement f b = (f x f y f z ) T is related to the velocity 

increment measurements Av b = (av^ Av y Av r ) T with f b = Av b jAt , equation 41 can be 
given as: 

W\t M ) = RlAv b (t k+l ) - (2ai + 02,)K< (fc ) A/ + g £ A / (42) 

15 

where Qf g and Q e ei are given in equations 5 and 6 respectively, while g £ is the gravity field 
vector and is defined as shown in equation 3. 

The first term on the right-hand side of equation 42 is the measured velocity increments 
20 transformed to the £ -frame. The second term is the Coriolis correction that compensates for 
the Earth rotation and the £ -frame change of orientation. The third tern is the gravity 
correction. 
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The velocity components at the current time epoch t k+l are computed using the modified 
Euler formula as follows: 



r%+\)-v%)+\kr'itk) + *r\tk + i)) (43) 

where v e = (v e V" v u ) are the velocity components as described in the £ -frame. 

6. Determine the altitude of the bottom hole assembly, h , (i.e. the true vertical depth TVD) 

using the vertical component of the velocity vector, V u . Since h = V u according to equation 
2, the TVD can be computed at time instant t k+l using the modified Euler formula as 
follows: 



*(/*+! ) = Kh ) + ^(v u (/* + , ) + V (t k ))A/ (44) 



Once the TVD is determined, the latitude q> , and the longitude X , can be computed. 
According to equation 2, <p and X can be computed at time instant ) using the modified 
Euler formula as follows: 



1 ) = <P<tk ) + T ^ !Lu At (45) 

2 M + n 



2 [N + h)cos(p 



It should be noted that the computation of the longitude, X , will suffer from singularities if 
this procedure is performed closer to the poles where the latitude cp becomes equal to n/2. 
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The six steps described above are repeated for each new set of measurements received from 
the two FOGs and the three accelerometers until 20° inclination. A block diagram describing the 
functionality of this procedure is shown on Figure 14. 



5 In an alternative embodiment of the invention, surveying may be accomplished with only the 
toroidal FOG and the three accelerometers until 20°inclination. A block diagram describing the 
functionality of this procedure is shown in Figure 14A. In this case, the toroidal FOG may 
monitor both azimuth and roll angles because of the near vertically of the BHA. As the pitch 
angle increases, the unity between azimuth and roll decreases and the ability of the toroidal FOG 
1 0 to monitor roll angle diminishes . 

Surveying the radical section after 20° inclination 

As shown earlier on Figure 13, change of axes orientation may take place at 20° 
15 inclination in one embodiment of the invention. The Y-axis is now oriented along the tool spin 
axis while Z ] -axis is oriented along the orthogonal direction. The X-axis remains collinear to the 
old X-axis. The initial values for the navigation parameters (position components, velocity 
components and attitude components) for surveying this part of the radical section of the well can 
be obtained from the last values determined just before changing the axes orientation. It should 
20 be reiterated that the location of the two FOGs and the accelerometers would not physically 
change due to the change of axes orientation. Once the procedure for the first part of the radical 
section of the well (before 20° inclination) indicates an inclination of 20°, the output of the 
toroidal FOG is interpreted as A0 y , since its sensitive axis is reoriented along Y l -axis. Similarly, 

the output of the normal FOG is interpreted as A0 Z since its sensitive axis is reoriented along the 
25 Z^axis (see Figure 12). In addition, the output of the accelerometer installed along the tool spin 
axis (Y-axis) is interpreted as Av^ instead of Av z while the output of the accelerometer 

installed along the orthogonal direction (Z^axis) is interpreted as Av z instead of &v y . The third 
accelerometer installed along the X-axis has its output interpreted as Av x the same as in the first 
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part of the drilling of the radical section of the well 



The navigation procedure for this part of the radical section of the well is exactly the 
same as the procedure described in the previous section except for three changes. The first 
5 change is the different interpretation of the measurements provided by the two FOGs and the 
three accelerometers. The second change is that the tool face angle becomes no longer equivalent 
to the azimuth, Instead, it becomes equivalent to the roll. The third change is that the inclination 
is no longer equivalent to the pitch (9 ) and is given as k\2 - 6 . The navigation procedure for 
surveying after 20° inclination is shown schematically in the block diagram shown as Figure 15. 

10 

Surveying the horizontal section of the well: 



Surveying the horizontal section of the well is different from the radical section due to the 
rotary mode of drilling during which the whole drill pipe starts to continuously rotate about its spin 

15 axis. Therefore, continuous surveying can no longer be performed and station-based surveying 
procedure should be prepared. The station-based surveying is based on stopping the drilling 
operation at certain survey stations. Since the surveying setup is completely stationary at each 
station, the two FOGs are affected by the Earth rotation rate and the three accelerometers are affected 
by the Earth gravity field. The arrangement of the surveying instruments inside the bottom hole 

20 assembly is shown on Figure 1 6. The relationship between the Earth rotation rate vector monitored 

by the FOGs and the same vector defined at the i -frame is given using the rotation matrix R\ as 
follows: 



6)„ 



cos ^ cos ^ - sin y/ sin 9 sin <f) sin y/ cos <j) + cos yf sin 9 sin <j> - cos 9 sin ^ 

-siny/xostf cos y/ cos 9 sin 9 

cos y/ sin <f> + sin y/ sin 9 cos </> sin y/ sin <j> ~ cos y/ sin 9 cos </> cos 9 cos <f) 
(47) 
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The angular velocity <» ;. mn , 

y a y ls monitored bv th* ^ • j , 
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(48) 

The component of the Earth , • 
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^=(cos^co S ^ cos ^ (s . n ^ 
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mpuiation, equation 49 can be 
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To determine the azimuth angle from equation 50, the pitch angle, 0 , should be 
computed. Moreover, complete station based surveying requires the determination of the 
inclination angle 0/2 - 0 ) and the tool face angle (the roll angle, § ). This information is 
determined by manipulating the accelerometer measurements. The accelerometer measurements 
{fx fy fz) T 3X0 related to the Earth gravity filed through the rotation matrix . This 
relationship can be written as: 





(f x ) 


f 


g b = 


fy 


-K £ g - 




\fzj 


V 



cosy cos 0 -siny sin#sin^ siny cos <f> + cosy sin 0 sin <f> -cos#sin^^ 
- sin y/ cos 9 cos y cos 9 sin 9 



(51) 



The expressions for f x ,f y and f z can be expressed in terms of the pitch and the roll angles by 
manipulating equation 5 1 as follows: 



f x = (cos 9 sin </)) y 
f y =-(sm0) y 
f z =-(cos9cos<p)y 



(52) 
(53) 
(54) 



From equations (52-54) the pitch and the roll angles can be expresses as: 



( I —\ 



9 = arctan 
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</> = arctan 



V fz J 



(55) 



(56) 
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The specific force measurements f x ,f y and f z are obtained from the velocity increment 
outputs of the accelerometers Av^, Av y and Av 2 as follows: 



10 



Av Av v Av 



Therefore, station-based surveying during the horizontal section of the well delivers only 
the orientation of the bottom hole assembly using equations 50, 55 and 56 to provide the 
azimuth, the pitch and the roll respectively. The inclination is then equal to njl - 0 and the tool 
face is exactly the roll. 

Error analysis for the station-based surveying: 



Although the measurements from the toroidal FOG and the accelerometers are compensated 
for all bias terms, there are still some residual errors that affect the accuracy of the attitude angles. 
1 5 These residual errors include the cross-axis sensitivity, the bias calibration error, the bias temperature 
shift, and the scale factor calibration error. Some typical values are assigned to these errors and they 
are defined by the manufacturer. 



The error components in the specific force measurements can be obtained by differentiating 
20 equations 52, 53 and 54. 

8f x = (-sin0sin^<5# + cosO cos <j> 8<fi)y (58) 
df y = - (cos 080) y (59) 
8f z = (sin 0 cos (j>80 + cos 0 sin^ 8(f) y (60) 

25 

The error in the pitch angle, 59 , can be obtained from equation 59 and is expressed as follows: 
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50 = - 



cos# 



(61) 



Using equations 58 and 60, {<5f x ) Z + (Sf z ) z can be given as follows: 



Wx? + Wzf = (sin 2 9 [S6f + cos 2 0 (<^) 2 )/ 2 



(62) 



Equation (62) can be manipulated with the substitution of the error in the pitch angle from 
equation (61) to give the error in the roll angle 8<|> as follows: 



10 



cos# 



(63) 



The errors in both the pitch and the roll are determined from equations 61 and 63 and are 
used to correct the values of the pitch and the roll angles obtained by equations 55 and 56. The 
correct value of the pitch angle is then used in equation 50 to determine the azimuth. However, 
1 5 any residual errors existing in the toroidal FOG affect the azimuth accuracy. The azimuth error 
can be determined from equation 49. If this equation is differentiated with the assumption that 
the pitch angle is free from error, the toroidal FOG error can be expressed as: 

8ca y - 5 (cosy/) cos 0{a> e cos cp) (64) 

20 

Then the error term 8 (cos\|/ ) can be written as 
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Sco, 



<?(cos^ l(ft/c ° S ^ 
cos# 



(65) 



As discussed above regarding the establishment of the desired azimuth direction, the error 
in the cosy term should be determined using equation 65 and completely removed before the 
value of the azimuth angle is computed. This is to avoid the magnification of the azimuth error if 
calculated directly as follows: 



Sy/ 



So, 



(co e cos^) y 
cos# 



smy/ 



(66) 



1 0 It can be inferred from equation 66 that the term 1/sin y/ largely magnifies the azimuth 

error, 8v|/ , as the bottom hole assembly approaches the North-South direction (i.e. when vf * 0 ). 
In addition, singularity problems exist in the calculation of 8v|r if v|/ = 0 . Therefore, to avoid all 
these problems, the error term 5 (cosy ) is computed prior to determining the azimuth. 



1 5 The whole station-based surveying procedure can be shown as a block diagram in Figure 

17. It should be mentioned that this procedure gives coarse computation of the attitude angles. 
Precise computation should be applied to improve the estimation accuracy. The technique for 
precise computation is described in the following sections. 



20 Error analysis of the motion of the bottom hole assembly 



In practical systems employing the technology of strapdown inertial navigation the 
accuracy is limited by the errors originating in the raw measurements as well as by the 
imperfections in the various components which are combined to build the system. These error 
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sources are categorized as inertial sensor errors and computational errors [Stephenson et al., 
1992; Schwarz et al, 1990]. Any lack of precision in the measurements of an inertial navigation 
system is passed from one estimate to the next with the overall uncertainty in the precision of the 
estimated quantity drifting with time. This explains why the accuracy of the INS deteriorates in 
5 the long term. For MWD surveying, the drilling process continues for several hours. Therefore, 
error models are required for analysis and estimation of the various error sources associated with 
the inertial system proposed in the previous section. 

Errors in dynamic systems are variable in time and therefore described by differential 
10 equations. Linearization of the non-linear dynamic system is the most common approach for 

deriving a set of linear differential equations which define the error states of the dynamic system. 
The error states expressed in the t -frame for the mechanization equations given in the previous 
section are defined as follows: 



15 e* =(5<pM,Sh 9 5V e9 SV n ,SV u ,S0,S<f>,Siff (67) 
where 

5cp ,5 A,, 5 A are the coordinate errors; 

SV e , 5V n , 5V u are the velocity errors ; 

59 , 5<|> , 5v|/ are the attitude errors; 

5co x , 8co y , 8co z are the errors in the angular velocity measurements; 

tifx&yi $fz are toe errors in the specific force measurements. 



Attitude errors 



20 The attitude errors are defined as the misalignment due to the orthogonal transformation 

between the b-frame and the I -frame. These errors (50 , 5<() and 5i|/ ) depend on two main 
sources. The first source are the errors in the measurement of angular velocities provided by the 
FOGs. The second source are the errors generated due to the Earth rotation and the change of 
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orientation of the i -frame. The expression describing the time rate of change of the attitude 
errors is given as follows: 

e l = -Q £ ie £ e - Sa)f e + R e b d (68) 

where 

e e ={S0 8<t> S\f/) T is the attitude error; 

= Qj + is the skew symmetric matrix for the rotation of the I -frame with respect to 

the i-frame; 

p p p p 
Sco i( = 5(o ie + So) e( is the angular velocity error vector corresponding to co u ; 

d = {dco x Sco y 8(0 z Y * s me vector of angular velocities measurement errors (FOG drift rates). 

The second term in the right hand side of equation 68, Sco^ , depends implicitly on the 
velocity errors SV e = (sv e SV n SV") T and the coordinate errors Sr l = (S<p SA 5h) T . 
This angular velocity error consists of two terms, Scof e and 5co ei , which are given as: 



8(o\ e = 



0 0 0^ 
- co e sin (p 0 0 
^ co e coscp 0 0^ 



8X 



(69) 
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(70) 
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With the substitution of the expressions of 8co\ e (equation 69) and Sco^ (equation 70), the 
attitude error equation (equation 68) can be rewritten as: 



co smcp + 



e . F tanffl 

co sin cp 

N + h 



co costp + 



N + h 



V'tancp 
N + h 
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-co cos<p- 
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M + h 
0 
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81 
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v J 



(71) 



Coordinate errors: 



10 Since the coordinates of the bottom hole assembly depend only on the velocity 

components V e , V n and V u , the coordinate errors Sep , bX mdbh are given directly in terms of 

velocity errors 5V e ,5 V n and 5V U . The expression for the coordinate errors is obtained by 
differentiating both sides of equation 2 that describes the relationship between the time rate of 
change of the coordinate components and the velocity components. Since this relationship is 

15 given by the matrix D~ ] with its elements depending on the coordinate components 9 and h , 
the coordinate errors 5cp, 6X mdbh are written as follows: 
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Sr ( = 
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(72) 



Velocity errors: 

The complete implementation of the £ -frame mechanization equations requires the 

computation of the three velocity components V e , V n and V u using equation 4. The 
corresponding error equations are obtained by differentiating equation 4 and are given as: 
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r 5V^ 
5V n 

sv" 

\ J 



where 
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(73) 



f e ,f„mdf u 

are three components of the specific force vector / -{f e f n f u ) T = Rbf 
= Rb(fx f y fzY 

, where f X9 f y and f z are the specific force measurement 



and R 



outputs of the accelerometers. 
is the mean radius of the Earth. 



Measurement errors: 
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The sensor measurements can be compensated for the constant biases in both the FOGs 
and the accelerometers. However, there are still some residual errors which are random in nature. 
These errors become significant in applications that require long-time surveying. Therefore, these 
residual random errors are modelled as stochastic processes. The residual gyroscopic 
5 measurement errors are correlated in time, therefore they are usually modelled as 1 st order Gauss- 
Markov processes as follows [Brown et al., 1992; Gelb, 1979]: 
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(74) 



where 

a x , a y and a 2 



are the reciprocals of the time correlation parameters of the random processes 
associated with the angular velocities measurements; 



a lx , <j ]y and <r ]z are the standard deviations of these random processes; 



W(t) 



is unity- variance white Gaussian noise. 



10 



Similarly, the residual accelerometer measurement errors are modelled as 1 st order Gauss- 
Markov processes as follows: 
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(75) 



15 where 

J3 x ,/3 y and f5 z are the reciprocals of the time correlation parameters of the random processes 

associated with the specific force measurements; 
<J 2x> a 2y m ^ <7 2z are ^ e stan dard deviations of these random processes. 
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Equations 71, 72 and 73 for the attitude, position and velocity errors are combined 
together with equations 74 and 75 for the FOGs and accelerometer errors to describe the behavior 

of the error state e i defined in equation 67. The augmented equation is given as a linear system 
with random system noise as follows: 

e e = Fe £ +GW(t) (76) 

where 

F is a 15 -by- 15 matrix known as the dynamic matrix of the error states; 

G is a 15-by-l vector known as the vector associated with the system noise; 

W{t) is unity-variance white Gaussian noise. 

The linear dynamic system shown in equation 75 can be modelled by the following 
discrete linear state equation: 

= F k^k e k + G k W k (77) 

where 

e k+l is the error state vector at time t k+] ; 
e k is the error state vector at time t k ; 

F k+lk is the transition matrix of the discrete linear system determined from the dynamic matrix of 
the error state F as F k + ] k = / + FAt , where / is the identity matrix and At = t k+l - t k ; 

W k is the random sequence of system noise. 

The problem now is to estimate the error state e k using the discrete linear state equation 
given in equation 76. Due to the random noise component existing in this linear state equation, 
its solution becomes a filtering problem which is usually addressed using discrete Kalman filter 
theory. 
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Aided strapdown inertia! navigation employing Kalman filtering theory: 

Kalman filtering is a well-known technique for applied optimal estimation especially in 
5 the field of inertial navigation and kinematic positioning [Brown et al., 1992]. Kalman filtering is 
applied to enhance the performance of the developed techniques by removing the effect of 
residual random errors during continuous surveying. In addition, Kalman filter delivers real-time 
statistical data related to the accuracy of the estimated values. To apply the Kalman filter, an 
independent measurement equation should be specified in addition to the system equation 
10 described in the previous section (equation 76). This measurement equation should be based on 
an independent measurement of superior accuracy different from the measurements provided by 
the FOGs and the accelerometers [Hulsing, 1989]. 

In MWD borehole surveying, continuous measurement of the length of the drill pipe as 
1 5 well as real-time measurement of the BHA penetration rate can respectively provide information 
about the depth of the borehole assembly (i.e. the variable h in the navigation vector) and the 
three velocity components V\V n and V u . This independent measurement provides what is 
known as aided strapdown inertial navigation. Figure 18 shows a block diagram of the technique 
of aided inertial navigation employing Kalman filtering. 

20 

The measurements provided by the FOGs and the accelerometers are sent to the surface 
of the drilling site with an appropriate telemetry technique. These measurements are processed 
with the surface computer system to provide the navigation parameters (three position 
components, three velocity components and three attitude angles). The error states are then 
25 estimated using Kalman filtering which utilizes the measurement of the length of the drill pipe 
and the BHA penetration rate as an independent information of superior accuracy. According to 
the technique of aided inertial navigation [Brown and Hwang, 1992], the length of the drill pipe, 

L , should be converted first into the corresponding true vertical depth (TVD), h* , as follows: 
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h* = L cos^ (78) 

where 

£ Is the inclination angle. 

/j* Is the true vertical depth as calculated using the length of the drill pipe. 



The true vertical depth determined using the length of the drill pipe, h* , is compared to 
5 the corresponding value obtained by the computational navigation procedure, h . Similarly, the 
velocity components determined using the measurements of the BHA penetration rate are 
compares with the corresponding values obtained by the computational navigation procedure. 
The error in between is processed by Kalman filter to provide an estimate of the error states e 
which is removed from the output of the navigation procedure to deliver corrected inertial output. 
10 Therefore, the measurement equation based on the length of the drill pipe as an aiding source is 
expressed mathematically as follows: 

y k =H k e k +v k (79) 

where 

y k is the error in the measurement vector of both the drill pipe length and the penetration rate. 
H k is the design vector of the measurement equation. 
v k is the discrete measurement random noise. 

15 

Since the design vector, H k , relates the measurement error in the TVD and the velocity 
components to the error vector, e k , it can be written as 
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The problem now is to estimate the error states e k based on the set of independent 
information y k with the assumption that the system noise W k and the measurement noise v k are 
uncorrected. The covariance matrix of the system noise Q k and the covariance matrix of the 
5 measurement noise R k are given as 



10 



20 



*riH'* III M 



with £[^v[]= 0 for all values of k and I . 



The Kalman filter provides a sequential recursive method for the optimal least mean 
variance estimation of the system states e k . The Kalman filtering is a two step process. The first 

step is to make a prediction of the error states based on the estimation of the previous error states. 
The second step is to update the predicted estimation using the available measurement from the 
15 current step. This approach can be quantified as [Schwarz et al., 1999]: 



(83) 



Prediction 

«*(-) = F kj c -\h-\(+) 

Where 

e k (-) is the predicted estimate of the error state at time t k ; 

e k _ x (+) is the updated estimate of the error state at time t k _ x ; 

/*(-) is the predicted estimate of the covariance matrix of the estimation error 

{p k {-) = E{{e k {-)-e k }[e k {-)-e k ] T )). 
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Update 

h (+) = h H + K k [y k - H k e k (-)] 

K k = P k {-)H T k \H k P k {-)H T k +R k Y (84) 
P k (+) = [l-K k H k ]/>,(-) 

where 

^ (+) is the updated estimate of the error state at time t k ; 

is the Kalman gain matrix; 
P k (+) is the updated estimate of the covariance matrix of the estimation error 

{p k (+) = E{[e k ( + )-e k ][e k (+)-e k ] T }) 

5 It should be noted that aided inertial navigation based on the length of the drill pipe as an 

aiding source will improve the accuracy of the navigation parameters on the vertical channel in a 
more direct way. The parameters of the vertical channel include the true vertical depth h , the 

vertical velocity V u . On the other hand, the accuracy of all other navigation parameters 
(horizontal channel parameters) will not be significantly affected by the length of the drill pipe as 

10 an aiding source. Alternatively, the independent measurement of the BHA penetration rate will 
improve the accuracy of the velocity components V e and V n and consequently the inclination 
and roll angles. In addition, the navigation parameters of the horizontal channels (East and North 
channels) have another advantage. The errors in the navigation parameters of the horizontal 
channels are bound due to the Schuller effect and oscillate with a period of 84 minutes (a 

15 frequency of about 1/5000 known as the Schuller frequency) [Mohammed, 1999; Salychev, 
1998]. 



As mentioned earlier, the technique of aided inertial navigation based on independent 
measurement of superior accuracy suits the continuous surveying procedure performed inside the 
20 radical section of the well. During the station-based surveying, the attitude angles are determined 
using the procedure shown on Figure 17. The values determined using this procedure are 
considered coarse computation and precise computation is still needed to improve the estimation 
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accuracy. The technique of precise computation is known as zero velocity update (ZUPT). The 
ZUPT technique utilizes the condition of zero velocity existing at each survey station (while the 
complete surveying setup is stationary) to update the estimate of errors in the Kalman filtering 
method. 

5 

Precise computation of the attitude angles during station-based surveying employing 
Kalman filtering theory: 

The augmented linear state equation (equation 76) describing the behavior of the error 
10 states in the £ -frame and its corresponding discrete representation is used for station-based 

precise computation. The only change introduced to suit the station-based surveying is to set the 
velocity components V e , V and V equal to zero. Therefore, the position error state equation is 
given as follows (see equation 72): 



15 
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(85) 



Similarly, the velocity error state equation and the attitude error state equation are given 
respectively as follows (see equation 73 and 71): 
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and 
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The measurement error states of both the accelerometer and the FOGs are described as 
first order Gauss-Markov process and they are expressed as follows: 
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(89) 



The augmented error linear state equation during the station-based surveying is obtained by 
1 5 combining equations (85-89). This augmented equation is similar to equation 77 but with 

different definitions of the dynamic matrix F and the system noise vector G which are given as: 
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and 

5 

G=(0 0 0 0 0 0 0 0 0 4 2a x°\x fay°\y &<*z<*\z ^Px°2x \Wy°\y ^JAz^ 

(91) 



In the case of station-based surveying, the measurement equation, y k = H k e k + v k , is 

10 completely different from the aided inertial navigation technique discussed in the previous 
section. There is no aiding source that can provide a measurement of superior accuracy. 
Alternatively, we employ the condition of zero velocity to construct the measurement equation. 
This condition is introduced because the bottom hole assembly is at a complete rest at each 
survey station. In these circumstances the error measurement vector y k is given as follows 

15 

y k =V e -V° (92) 

where 
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V 1 is the velocity vector expressed in the t -frame, V e =\V e V n V u 
V° is the zero velocity vector, V° = (0 0 0) 7 . 



Since the error measurement vector, y k , expresses an error in velocity, the measurement 
equation is given similarly to equation 79 but with the following definition of the design matrix 



H = 



fo 0010000000000 0^ 
000010000000000 
000001000000000 



(93) 



Kalman filtering described in equations 83 and 84 can then be applied but using the new 
definitions for F and H (the dynamic and the design matrices respectively), as well as the new 
10 definitions of the G and y k (the system noise and error measurement vectors). 



The precise computation technique described in this section improves the estimation 
accuracy by removing the residual random noise components existing in the accelerometer and 
the FOG measurements. Since the ZUPT technique is directly related to the velocity components 

15 of the bottom hole assembly, the velocity errors, SV e ,S V n and 5V U , will be significantly 
affected in a more direct way than the position errors, 5q> , h\ and hh , and the attitude errors, 
80 ,5<|) andSij/ . The precise computation usually takes about 10 to 15 minutes for the Kalman 
filter parameters to converge. The exact time needed for convergence depends mainly on the 
sensor accuracy and the level of the random noise at their outputs. 

20 

As will be apparent to those skilled in the art, various modifications, adaptations and 
variations of the foregoing specific disclosure can be made without departing from the scope of the 
invention claimed herein. 
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